Transverse Velocity Dependence of the Proton- Antiproton Ratio as a Signature of the 

QCD Critical Point 
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The presence of a critical point in the QCD phase diagram can deform the trajectories describing 
the evolution of the expanding fireball in the hb — T phase diagram. If the average emission time 
of hadrons is a function of transverse velocity, as microscopic simulations of the hadronic freeze- 
out dynamics suggest, the deformation of the hydrodynamic trajectories will change the transverse 
velocity (/Jr) dependence of the proton-antiproton ratio when the fireball passes in the vicinity of 
the critical point. An unusual /^-dependence of the p/p ratio in a narrow beam energy window 
would thus signal the presence of the critical point. 



Lattice-QCD simulations have shown that the tran- 
sition between the hadronic and quark-gluon plasma 
phases of quantum chromodynamics (QCD) at vanishing 
baryon chemical potential \xb is a crossover transition 
This raises the question whether the crossover transition 
becomes a first-order phase transition for larger values 
of [ib- Several attempts have been made to locate the 
critical point, i. e. the endpoint of the first-order transi- 
tion line, in lattice simulations 0,0, 0], but its existence 
is still in doubt [8]. The inconclusive theoretical results 
have motivated plans for a systematic exploration of the 
properties of hot QCD matter as a function of the net 
baryon density by means of a collision energy scan at the 
Relativistic Heavy Ion Collider (RHIC) [fj|,|6|. The search 
for the QCD critical point also forms part of the motiva- 
tion for the NA61 experiment at the CERN-SPS and 
for a new facility dedicated to the study of compressed 
baryonic matter at the Facility for Antiproton and Ion 
Research (FAIR) in Germany. 

Ideas for experimental signatures for the presence of 
the critical point have mostly focused on fluctuations in 
certain observables related to the order parameter of the 
chiral transition 0, [l(| • General arguments lead one to 
believe that such fluctuations are enhanced in the vicin- 
ity of the critical point. Unfortunately, several reasons 
throw doubt on the usefulness of fluctuation observables 
as practical signatures of the QCD critical point. First, 
fluctuations are suppressed, compared to the static case, 
when the matter passes rapidly through the critical re- 
gion during the expansion due to critical slowing down 
[111 ]. Secondly, the hot matter does not freeze out at the 
critical point, but at a much lower temperature, when 
the critical fluctuations may well have been washed out. 
Finally, it is unclear in which observable fluctuations are 
most promising experimentally. In exploratory experi- 
ments at the CERN-SPS, only fluctuations in the K/tt 
ratio at beam energies below 40 GeV/A have shown signs 
of an unusual behavior [T2, EH ■ 

Here we propose a possible signature of the presence of 



a critical point in the QCD phase diagram, which may be 
more robust than fluctuations associated with the order 
parameter of the chiral phase transition. Our idea is 
based on the observation that the critical point serves 
as an attractor of the hydrodynamical trajectories in the 
fis — T plane describing the expansion of the hot matter 
|14j . We describe below how this focusing effect manifests 
itself in an experimental observable. 

The universality argument tells us that the critical ex- 
ponents around second order phase transitions are de- 
termined only by the dimensionality and symmetry of 
the system. The QCD critial point, if it exists, belongs 
to the same universality class as the 3-dimensional Ising 
model and liquid-gas phase transition [l(|. The singu- 
lar part of the thermodynamic variables near the crit- 
ical point is a function of two variables, which can be 
mapped onto the variables characterizing the phase dia- 
gram of the 3-dimensional Ising model: the reduced tem- 
perature r = (T — T c ) /T c and the external magnetic field 
h. In the QCD phase diagram, the axis corresponding 
to r points in the direction of the phase boundary; the 
direction of the axis associated with the variable h is not 
known ll|, |l5|. However, it is clear that the critical re- 



gion is more elongated along the r-direction, because the 
critical exponent associated with r is larger than that 
associated with h 



11] 



The focusing effect can now be understood as follows. 
The entropy density s and the baryon density depend 
in different ways on r and h. As a result, the ratio s/iif,, 
which is constant along an isentropic trajectory, assumes 
many different values in the vicinity of the critical point. 
Therefore, hydrodynamic trajectories for a range of dif- 
ferent values of s/ni, pass near the critical point, thus 
causing the focusing effect. 

The extent of the focusing region depends on the size 
of the critical region in the /j^ — T plane, in which ther- 
modynamic susceptibilities are significantly enhanced by 
the critical exponents. The size of the attractive basin 
can, in principle, be determined by lattice-QCD Simula- 
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tions. At the moment this information is not available, 
as the location and even the existence of the critical point 
in QCD are not established. Model studies in simplified 
theories suggest that the size of the attractive region is 
sensitive to calculational details jlq ]. Here we will sim- 
ply assume that the critical region is sufficiently large to 
induce a significant focusing effect. We use the model of 
Nonaka and Asakawa [3] to describe the influence of the 
critical point on the thermodynamic variables. In this 
model the entropy density is obtained by interpolation 
between the entropy densities sh, sq of the hadronic and 
quark phase as 

s(T,/j, B ) = i (1 - tanhS'c) s H (r, Mb) 

+ i(l + tanhS c ) SQ (l> B ), (1) 

where S c (T,iib) is proportional to the critical part of 
the entropy density obtained by rescaling the expressions 
from the 3-dimensional Ising model. The proportion- 
ality constant entering into the definition of S c deter- 
mines the size of the influence region of the critical point. 
Here we differ from ref. [14j by choosing the parameters 
AT crit = 20 MeV,A/i cr it = 100 MeV,L> = 0.5. This 
choice yields a much narrower critical region as shown 
in Fig. Q] The width of the cross-over between the two 
phases at /is = is approximately 45 MeV, in rough 
agreement with lattice-QCD results 17 1. 
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FIG. 1: Critical region in the fiB — T plane. The thick 
solid line to the right of the black dot shows the location of 
the phase boundary, starting at the critical point (fj, c ,T c ) — 
(550, 159) MeV. The thin lines to the left indicate the contours 
of equal values of the cross-over parameter tanh S c between 
—0.8 and +0.8 in increments of 0.2. 

The main characteristic of the fireball evolution in the 
presence of a critical point is that hydrodynamical tra- 
jectories, which would normally tilt to the right after 
crossing the phase boundary (see solid line in Fig. fj] for 
a smooth cross-over or the dash-dotted line for a first- 
order transition) , make a detour into the vicinity of the 



critical point and then turn to the left as the tempera- 
ture falls below T c (see dashed line in Fig. [2]). For our 
argument, the important difference is the behavior just 
below T c , where both T and /j,b decrease for the trajec- 
tory through the critical point, while /is stays roughly 
constant or increases with falling temperature for trajec- 
tories away from the critical point. This difference can 
have visible consequences if hadron emission occurs over 
a finite range of temperatures, and if emission from differ- 
ent points along the trajectory can be discriminated. For 
instance, the ratio hb/T monotonically increases below 
T c along the "normal" (solid or dash-dotted) trajectories 
in Fig. [21 implying a falling antiproton-to-proton (p/p) 
ratio. On the other hand, the dashed trajectory in Fig. [2] 
implies an approximately constant or even slightly de- 
creasing value of (J-b/T and thus a rising p/p ratio as the 
temperature falls below T c . 
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FIG. 2: Hydrodynamical trajectories in the QCD phase dia- 
gram with and without the presence of a critical point. Possi- 
ble trajectories in the /is — T plane in the absence of a critical 
point are shown as solid line (for a cross-over transition) or 
dash-dotted line (for a first-order transition); the trajectory 
in the presence of a critical point is shown as dashed line. 
All trajectories meet at the bulk chemical freeze-out point. 
Arrows indicate the direction of time evolution. 

In order to confirm this qualitative argument we 
present a quantitative analysis based on the assumption 
that the attractive basin of the critical point is reached in 
central Pb+Pb collisions at 40 GcV/A. In Fig.[3]we show 
the p/p ratio along the three trajectories shown in Fig. [5] 
as a function of the entropy density between T c and the 
chemical freeze-out point, which has been determined to 
lie at (^ ch , X^) « (400, 145) MeV by a statistical model 
fit to experimental data [22j |. As anticipated, the p/p 
ratio falls or remains constant between the phase bound- 
ary and chemical freeze-out for the "normal" trajectories 
(solid and dash-dotted lines), but rises for the trajectory 
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deformed by the presence of the critical point (dashed 
line). 
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FIG. 3: Antiproton-to-proton ratio along the trajectories 
shown in Fig. [2] as a function of the entropy density. The 
curves start at the phase boundary T c « 160 MeV and con- 
tinue down to T « 110 MeV. The location of the chemical 
freeze-out point fjUch,Tch) deduced from experimental data is 
indicated by the open and solid squares. Note that the p/p 
ratio fonly rises for the trajectory deformed by the critical 
point. 

We next discuss how baryon emission from different 
points along the hydrodynamical trajectory may be dis- 
tinguished. We first note that data from Au+Au colli- 
sions at RHIC have been explained by the assumption 
that the emission of hadrons with intermediate trans- 
verse momentum (pr ~ 2 — 5 GeV/c) occurs at the 
phase boundary by recombination of constituent quarks 
[18l [l9j . Bulk freeze-out of hadrons, on the other hand, 
occurs when the mean free path of hadrons becomes com- 
parable to the size of the fireball. The mean free path 
relevant to transport properties generally grows with in- 
creasing hadron momentum. This implies that hadrons 
with large transverse momentum should freeze out ear- 
lier, on average, than hadrons with a small transverse mo- 
mentum. In the extreme, intermediate pt hadrons may 
be produced at or near the phase boundary. This effect 
can also be understood by invoking detailed balance. A 
highly energetic hadron, impinging onto the fireball from 
the outside, would penetrate deeper into the matter than 
a low-energy hadron. Conversely, energetic hadrons will 
be emitted, on average from deeper inside the matter and 
thus earlier than low-energy hadrons. 

The differential emission of baryons as a function of 
transverse momentum can be analyzed quantitatively 
in the framework of a microscopic hadron transport 
model, e. g. UrQMD [2(J 21 1. Such transport models 



nary hadronic reactions are commonly used to describe 
the freeze-out and break-up of the fireball produced in 
relativistic heavy- ion collisions into hadrons. Utilizing 
the UrQMD model, we have calculated central Au+Au 
collisions at a fixed target energy of 40 GcV/A, which 
may lead to conditions for which the matter passes near 
the QCD critical point. We then determined the last 
time of interaction in the medium ("emission time") for 
all final-state (anti-)protons. In order to discriminate 
between slow and fast hadrons, we choose two trans- 
verse velocity windows: Pt < 0.25 for slow hadrons and 
0.5 < (3t < 0.75 for fast hadrons, respectively. 
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FIG. 4: Top: UrQMD predictions for the emission time (last 
interaction time) distribution of protons and antiprotons in 
central Au+Au collisions at 40 GeV/ A The emission time 
distributions are shown separately for the kinematic windows 
/3t < 0.25 and 0.5 < j3r < 0.75. Fast particles are emitted on 
average significantly earlier than slow particles. Bottom: p/p 
ratio as a function of transverse velocity fir. 



The distribution of emission times (counting from the 
moment of full overlap of the colliding nuclei) for fast 
(anti-)protons is shown in the top frame of Fig. [4] com- 
pared to the emission time distribution for slow (anti- 
protons. The emission of protons and antiprotons in the 
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0.5 < j3r < 0.75 window is clearly correlated with early 
emission times and occurs approximately 4-5 fm/c ear- 
lier than for slow (anti-)protons. The average emission 
times for fast protons and antiprotons are 19.5 ±0.02 and 
16.5±0.06 fm/c, respectively compared to 22.9±0.01 and 
21.2 ±0.2 fm/c for slow moving protons and antiprotons. 

The bottom frame of Fig. [4] confirms the expectation 
from the solid line in Fig. [3j that the p/p ratio should 
rise as a function of (3t in the absence of a critical point. 
Because the UrQMD calculation does not include finite- 
density corrections to the p+p annihilation cross section, 
the overall value of the ratio should not be compared with 
experiment. We also note that UrQMD does not contain 
any physics related to the QCD phase transition and only 
serves here as a model to study the correlation between 
emission time and transverse velocity. We expect the ob- 
served correlation to persist in the presence of a critical 
point. The fall in the p/p ratio for (3t > 0.8 is due to kine- 
matical constraints on the p production and independent 
of collision size and energy. Since only a very small frac- 
tion of the baryon yield resides in the (3t > 0.75 range we 
have selected the 0.5 < /3t < 0.75 range as representative 
for high /?t and early freeze-out nucleons. 

A search of existing data revealed that the p/p ratio in 
Pb±Pb collisions at the CERN-SPS has been measured 
as a function of beam energy by the NA49 collabora- 
tion, which has published proton and antiproton spectra 
for fixed-target beam energies of 20, 30, 40, 80, and 158 
GeV/A [23j ]. Interestingly, the antiproton spectrum mea- 
sured at 40 GeV/^4 exhibits an anomaly. Whereas the ex- 
ponential slope of the antiproton spectrum is flatter than 
the slope of the proton spectrum at other beam energies, 
it is slightly steeper at 40 GeV/A. A flatter antiproton 
spectrum is compatible with differential chemical freeze- 
out on a trajectory similar to the solid or dash-dotted 
trajectories shown in Fig. [2 while a steeper spectrum 
would require a trajectory of the type expected in the 
vicinity of the critical point (similar to the dashed line in 
Fig. [2]). The relative suppression of antiprotons at large 
transverse momentum is clearly visible in the spectrum 
itself (see Fig. 3 in ref. 23]). The size of the statistical 
errors of the measurement does not permit a firm con- 
clusion about this anomaly, but it is certainly compatible 
with the arguments presented here. 

In summary, we have shown that the evolution of 
the p/p ratio along isentropic curves between the phase 
boundary in the QCD phase diagram and the chemical 
freeze-out point is strongly dependent on the presence or 
absence of a critical point. When a critical point exists, 
the isentropic trajectory approximately corresponding to 
hydrodynamical expansion is deformed, and the p/p ra- 
tio grows during the approach to chemical freeze-out. If 
nucleons of high transverse momentum are chemically 
frozen out earlier than the slow nucleons, as it is sug- 
gested by microscopic simulations of hadronic dynamics, 
this result will translate into a p/p ratio that falls with 



increasing transverse momentum instead of a rise or flat 
behavior in scenarios without critical point. This behav- 
ior would only occur at those beam energies, for which 
the fireball reaches the critical point. Depending on the 
actual size of the attractive region around the critical 
point, the search for an anomaly in the j/t dependence of 
the p/p ratio may require small beam energy steps. Note 
that the location of the critical point in our model study 
was chosen such that it is encountered by the hydrody- 
namical trajectory for conditions reached for a beam en- 
ergy of 40 GeV/A and a fixed-target. For a different lo- 
cation of the critical point, similar behavior would occur 
at other beam energies. 
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